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ABSTRACT 



The meridional flow in the Sun is an axisymmetric flow that is generally poleward directed at 
the surface, and is presumed to be of fundamental importance in the generation and transport 
of magnetic fields. Its true shape and strength, however, is debated. We present a numerical 
simulation of helioseismic wave propagation in the whole solar interior in the presence of a 
prescribed, stationary, single-cell, deep meridional circulation serving as a test-bed for helioseismic 
measurement techniques. A deep-focusing time-distance helioseismology technique is applied to 
the artificial data showing that it can in fact be used to measure the effects of the meridional flow 
very deep in the solar convection zone. It is shown that the ray-approximation which is commonly 
used for interpretation of helioseismology measurements remains a reasonable approximation even 
for the very long distances between 12° and 42° corresponding to depths between 52 and 195 Mm 
considered here. From the measurement noise we extrapolate that on the order of a full solar 
cycle may be needed to probe the flow all the way to the base of the convection zone. 

Subject headings: Methods: numerical, Sun: helioseismology, Sun: interior, Sun: oscillations 



1. Motivation and Objectives 

The meridional circulation is known as a pole- 
ward flow near the solar surface and is believed 
to be an important part of the dynamo process in 
the solar convection zone. Its true strength and 
shape in the deeper interior is presently unknown 
or at least highly debated. There are conflict- 
ing theories and observational evidences. On the 
one hand, theoretical works (e.g., Kitchatinov & 
Riidigcr 2005) usually favor deep meridional flows. 
Due to mass conservation, the return flow then is 
relatively strong with flow velocities of the similar 
same order as the flow near the surface. 

On the other hand, observational evidence 
points towards a shallow meridional flow. Mitra- 
Kraev & Thompson (2007) measured p-mode fre- 
quency shifts from three months of SOHO/MDI 



data and inferred a meridional flow that reverses 
at a depth of about 40 Mm. They also find evi- 
dence of a possible second reversal deeper below. 
Hathaway (2011) uses the advection of convection 
cells by the meridional circulation to measure the 
flow velocity at depth and finds a return flow start- 
ing at a depth of only 35 Mm. Other techniques, 
such as time-distance helioseismology have also 
been used to measure the meridional flow (e.g., 
Giles et al. 1998; Chou & Dai 2001; Zhao & Koso- 
vichev 2004; Zhao et al. 2012) but primarily for 
small depths. 

The goal of this work was to numerically sim- 
ulate the propagation of helioseismic waves in the 
Sun in the presence of meridional and other large 
scale flows, and to use the artificial data from such 
simulations to evaluate the possibility of measur- 
ing small flows very deep in the Sun, in particular 
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using time-distance helioseismology. There have 
been previous works that addressed the measur- 
ability of deep meridional flows in the Sun such 
as Braun & Birch (2008). They concluded that 
as much observations as a full solar cycle may be 
needed to resolved flows near the base of the so- 
lar convection zone using the helioseismic holog- 
raphy technique. Various models of meridional 
flows have been proposed, and for this paper we 
have performed a numerical simulation of the solar 
acoustic wave field for a global-Sun model which 
included a meridional circulation with a deep re- 
turn flow, and have performed time-distance he- 
lioseismology measurements for this model as if 
it were observational data. We compare the mea- 
surements with the predictions of the ray-path he- 
lioseismology theory and estimate the noise level 
due to the stochastic nature of solar oscillations, 
and the sensitivity of the helioseismology tech- 
nique. 

2. Numerical Simulation 
2.1. Simulation Code 

We have built a numerical code that solves 
the linearized propagation of helioseismic waves 
throughout the entire solar interior in the presence 
of a background structure and flow model (Hartlep 
& Mansour 2005). This code has been used in pre- 
vious studies to simulate the effects of localized 
sound speed perturbations, e.g., for testing helio- 
seismic far-side imaging by simulating the effects 
of model sunspots on the acoustic field (Hartlep 
ct al. 2008; Ilonidis et al. 2009), for validat- 
ing time - distance helioseismic measurements of 
tachocline perturbations (Zhao et al. 2009), and 
for studying the effects of localized subsurface per- 
turbations (Hartlep et al. 2011). For the present 
case, the code has been extended to include the 
effects that mass flows have on the propagation of 
helioseismic waves. 

The simulation code models solar acoustic oscil- 
lations in a spherical domain using the Euler equa- 
tions linearized around a stationary background 
state characterized by the background density, p , 
mass flows velocity, v , sound speed, c , and accel- 
eration due to gravity, go. The equations derived 
for the perturbations around the base state are 



then: 

d t p' = -V-m' + S, (1) 

d t ra! = -Vc^' + p'go+voK,- Vp') (2) 

+ p'(v • Vv + v V • v ) - (v • Vm' 

+ m' • Vv + m'V • v + v V • m'), 

where p' and m' = p'v Q + p Q v' are the density 
and momentum perturbations associated with the 
waves, respectively. 

Several simplifications where used in deriving 
these equations. In particular, perturbations of 
the gravitational potential have been neglected, 
and the adiabatic approximation has been used. 
The entropy gradient of the background model 
has been neglected in order to make the linearized 
equations convectively stable. Previous calcula- 
tions have shown that this assumption does not 
significantly change the propagation properties of 
acoustic waves including their frequencies, except 
for the acoustic cut-off frequency, which is slightly 
reduced. Because of this simplification, no sepa- 
rate energy equation needs to be solved. 

The equations are formulated in a non-rotating 
frame. Rotation can however be accounted for by 
prescribing an appropriate flow. This approach 
saves computing the usual Coriolis and centrifu- 
gal forces that appear in the equations written 
in a rotating frame. There is no additional time- 
stepping or stability constraint on the numerical 
method because the rotation speed in the solar in- 
terior is always significantly smaller than the speed 
of sound. 

In the Sun, vigorous turbulent flows near the 
photosphere are the primary sources of acoustic 
perturbations. This however is a non-linear pro- 
cess and is lost by linearizing the equations. A 
random function S has therefore been added to 
Eqn. 1 to mimic the excitation of these acoustic 
perturbations. S is random in time and horizon- 
tal space, and peaks at a depth of 150 km below 
the photosphere. 

The equations are solved in spherical coor- 
dinates using a pseudo-spectral method. 2/3- 
dealiasing is used. Scalar quantities such as pres- 
sure and density are expanded in spherical har- 
monic basis functions for their angular structure 
and B-splines (de Boor 1987; Loulou et al. 1997; 
Hartlep & Mansour 2004: Hartlep & Mansour 
2005) for their radial dependence. Vector fields 
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such as the perturbation of the momentum are 
expanded in vector spherical harmonics and B- 
splines. Vector spherical harmonics (Hansen 1935; 
Hill 1954) arc analogous to spherical harmonic 
functions but for vector quantities. They were se- 
lected here because the coordinate singularities in 
spherical coordinates can be treated straightfor- 
wardly in these basis functions. Due to the na- 
ture of spherical harmonics, no special treatment 
is necessary for the poles, and only the center of 
the sphere needs special treatment. For the ex- 
panded variables to be finite at the center, we must 
enforce that the expansion coefficients follow spe- 
cific asymptotic behaviors as r — > 0. For details 
on how this is implemented using B-splincs, the 
reader is referred to Hartlep et al. (2006). 

B-splines of polynomial order 4 are used in our 
simulations. The spacing of the generating knot 
points is chosen to be proportional to the local 
speed of sound. This results in a high radial res- 
olution near the solar surface, where the sound 
speed is low (less than 7 km/s), and a signifi- 
cantly lower resolution in the deep interior, where 
the sound speed surpasses 500 km/s. This pro- 
vides a constant Courant-Friedrichs-Lewy (CFL) 
condition throughout the domain. A total of 350 
B-splines are used to discretize the entire simula- 
tion domain reaching from the solar center to an 
outer radius of 700 Mm. 

Effectively non-reflecting boundary conditions 
are used at the upper boundary by means of an 
absorbing buffer layer by adding terms —up' and 
— am' in the above equations. The equations are 
recast using an integrating factor expert, and ad- 
vanced using a staggered leapfrog scheme in which 
p' and m' are offset by half of a time step. The 
purpose of the buffer layer is to damp waves pass- 
ing thought the temperature minimum before they 
reach the numerical boundary. Such waves would 
ordinarily escape into the chromosphere, and we 
do not want them to artificially reflect back from 
the numerical boundary. The damping coefficient 
a is non-negative, smooth and constant in time. It 
is zero in the interior and increases smoothly into 
the buffer layer. Similar damping is used in the 
deepest interior near the solar center. Waves of 
high spherical harmonic degree do not travel very 
deep. Their lower turning radius r t given by 

cl(r t )/r^u 2 i '1(1 + 1), (3) 



and decreases with increasing I. to in the equation 
is the wave's temporal frequency. Therefore, as 
we go towards the center of the Sun, we only need 
to resolve waves of lesses and lesses spherical har- 
monic degrees. Carrying higher degrees than re- 
quired would, in fact, unnecessarily limit the time 
step. In order to avoid this, we use damping co- 
efficients that are /-dependent. Perturbations of 
spherical order I do not travel below their corre- 
sponding turning radius as given by Eqn. 3, and 
can be damped. However, we leave about 100 Mm 
of space below the turning radius before damping 
starts since the waves have finite extend. Such 
damping effectively removes high-/ perturbations 
in the deep interior without effecting the propaga- 
tion of helioseismic waves. Modes below I — 40 are 
not damped at all because wave of such / can ac- 
tually travel very close and even through the solar 
center. 

Without flows in the background model, waves 
can, within the physics of this setup, only gain en- 
ergy through the acoustic source term, and they 
can only loose energy through the top bound- 
ary. There is no feedback between waves and the 
sources. These simulations run stable and will 
eventually reach an energy equilibrium. In the the 
presence of mass flows in the base state, in particu- 
lar with strong velocity gradients, it is conceivable 
that waves may gain momentum and energy from 
wave-flow instabilities, such as Kelvin-Hclmholtz 
instability. Such growth would be proportional to 
the wave's own momentum. Our simulations have 
shown however that for solar-type flows such in- 
stabilities grow very slowly, and adding a small 
amount of viscous damping (not shown in the 
above equations) is enough to ensure long-term 
stability. 

2.2. Flow Model and Simulation Run 

For this study, we have simulated the propa- 
gation of helioseismic waves through a station- 
ary model of the solar meridional circulation. 
The original flow model is the "reference model" 
from Rempel (2006). A visualization of the model 
is shown in figure 1. This meridional flow is char- 
acterized by a single-cell circulation in each hemi- 
sphere with a deep return flow that starts at a 
depth of about 150 Mm below the surface. The 
model has a maximum poleward velocity near the 
surface of approximately 14 m/s and a maximum 
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Fig. 1. — Visualization of the meridional flow model 
of Rcmpcl (2006) used in the simulation. The arrows indi- 
cate the direction of the flow, and their size is proportional 
to the flow speed. The Sun's rotation axis is at the left 
edge of the panel with the north pole at the top. The flow 
in the southern hemisphere (not shown) is given by mirror- 
symmetry about the equatorial plane. 

return flow velocity of 3 m/s. Detecting such weak 
flows from helioseismic measurements is extremely 
challenging due to the inherent randomness of so- 
lar oscillations resulting in small signal-to-noise 
ratios (S/N). It is expected that very long obser- 
vations on the order of many years are needed 
to measure the flow in the deepest parts of the 
convection zone (Braun & Birch 2008). Unfortu- 
nately, it is not practical to simulate such long 
time series giving current computer capabilities. 
The present simulation required approximately 1 
day of computing time for every three hours of 
solar time using 264 cores on the Pleiades super- 
computer at the National Aeronautics and Space 
Administration's (NASA) Advanced Supercom- 
puting Division (NAS). However, the most mea- 
sured quantities in helioseismology - in particu- 
lar apparent travel-times, frequency shifts, etc - 
depend linearly on the causing perturbations in 
the solar interior, i.e. sound speed variations or 
flows. This is a true as long as the perturba- 
tions are small, i.e. in the case of flows: as long 
as the flow velocities are small compared to the 
speed of sound. Therefore, we can improve the 



S/N without changing the physics by uniformly 
increasing the amplitude of the background flow 
model in the simulation. Such a simulation run 
for a manageable time is in some sense equivalent 
to a simulation of the original model but run much 
longer - thanks to T.L. Duvall, Jr. (2011, private 
communication) for this suggestion. In this work, 
we have increased the meridional flow by a factor 
of approximately 36 such that the maximum flow 
velocity is 500 m/s. This speed is still significantly 
smaller than the sound speed anywhere in the so- 
lar interior, but provides a large increase in S/N. 
In fact, since the (uncorrelated) realization noise 
reduces approximately with the square root of the 
length of the time series, our simulation should 
have a S/N that is equivalent to a simulation of 
the original model but run for a period 36 2 times 
longer. 

In the following, we list other parameters of the 
simulation in the study. The simulation resolves 
spherical harmonic degrees from to 170. It is suf- 
ficient to capture only this range of small to mod- 
erate spherical harmonic degrees because we are 
only interested in large scale and deep flows. The 
time step was 1 second; results were saved with 
a cadence of 30 seconds; and the simulation pro- 
duced a total of approximately 76 hours of data. 
The first 4 hours of data where discarded because 
they represent transient behavior as the simula- 
tion was started from a model without any waves 
in it. In the end, a total of 4,096 minutes of data 
have been produced and analyzed. 

3. Measurement technique 

In this study, we want to measure the effects 
that the meridional flow in the simulated Sun has 
on acoustic waves. When waves are advected by 
a flow, their travel time is reduced when they 
travel in the same direction as the flow and their 
travel time increases when they travel against the 
flow. The travel times are also effected by local- 
ized sound speed variations, however, the differ- 
ence between the travel times of waves going in 
opposite directions along the same travel path is 
to first order not sensitive to local sound speed 
variations because both directions are effected the 
same way by such perturbations. The difference is 
only sensitive to the flow along the travel path. 

In the following we describe the scheme we used 
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Fig. 2. — Visualization of the deep focusing measurement 
scheme used in this study. Panel (a) shows examples of 
acoustic ray paths {solid lines) originating from and trav- 
eling to a range of annuli centered around a latitude of 30°. 
Panel (b) shows a measurement annulus and its decompo- 
sition into sectors. 

to measure the apparent acoustic travel-time dif- 
ferences between waves traveling northward and 
waves traveling southward, in the following re- 
ferred to as N-S travel-time difference, as well as 
the travel-time difference between east- and west- 
traveling waves (E-W travel-time difference) . This 
N-S difference is sensitive to the meridional flow. 

From the simulation, we used a 4,096-minute- 
long time series of the radial component of the 
perturbation velocity at a fixed geometric height 
of 300 km above the photosphere. The simula- 
tion provides data on the full solar surface, how- 
ever, we processed these data more similar to 
how one would process observations. First, the 
time series is split in 1,024-minutc-long segments 
with 50% overlap. In each segment, we select 
120° x 120° tiles centered at 0°, 90°, 180°, and 
270° in longitude and -30°, 0°, and +30° in lat- 
itude, and remap these into heliographic coordi- 
nates using Postel's projection with a pixel size of 
0.6°, the same resolution as SOHO /MDI mcdium- 
1 data (Kosovichev et al. 1997). 

The geometry of the measurement scheme is the 
following: On the surface, centered around each 
longitude and latitude, we select a series of one- 
pixel-wide annuli ranging in diameter from 12° to 
42° in 1.2° increments. Examples of ray paths 
traveling from a point on an annulus to an oppos- 



ing point is shown in Figure 2a. In each annu- 
lus, the instantaneous signal is averaged over 30°- 
wide sectors in the north, south, east, and west 
direction (Fig. 2b). We then cross-correlate the 
sector-averaged signal with that from its opposing 
sector, i.e. the north sector with the south sector 
and the west sector with the east sector. Cross- 
correlations for the same latitude within each tile 
are averaged together. The longitude-averaged 
cross-correlation functions for positive and nega- 
tive time-lag are then fitted separately using Ga- 
bor wavelets (Kosovichev & Duvall 1997), and the 
resulting phase-travel times are subtracted from 
each other yielding the apparent travel-time differ- 
ence between north- and south-going (N-S) waves 
and between east- and west-going (E-W) waves. A 
very small fraction of far outliers with travel-time 
differences larger than 2 minutes are misfittings, 
and are discarded. Measurement schemes of this 
type have previously been called deep-focusing 
schemes (e.g., Zhao ct al. 2009) but should not be 
confused with other deep-focusing schemes such 
as Duvall (1995) which are constructed such that 
ray paths cross or "focus" at a target location be- 
low the surface. Zhao et al. (2009) used a scheme 
very similar to the present one to measure deep 
sound speed perturbations. However, they em- 
ployed 90°-wide measurements sectors, i.e. quad- 
rants, instead of the 30°-wide sectors used here. 

4. Results 

4.1. Travel-time differences 

N-S travel-time differences measured from the 
simulation data for a range of distances are pre- 
sented in Figure 3. The travel distances of 12°, 
18°, 24°, 30°, 36°, and 42° correspond to lower 
turning points of the ray paths of approximately 
52, 78, 106, 134, 164, and 195 Mm below the pho- 
tosphere, respectively. Since individual measure- 
ments are fairly noisy, it is advantageous to trade 
some spatial resolution for reduced noise. The re- 
sults have therefore been averaged over a small 
range of travel distances (±2.4°) as well as a range 
of latitudes (±3°). 

Figure 3 also shows the N-S travel-time dif- 
ferences expected from ray-approximation calcu- 
lations for the prescribed flow. Assuming small 
perturbation, the ray paths are the same as in the 
case without flows, and the travel-time difference 
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Fig. 3. — Travel-time differences between north-going and south-going waves for 6 different travel distances of 12, 18, 24, 30, 
36, and 42 hcliographic degrees. The panels show the values measured from 4,096 minutes of simulation data (solid curves) and 
corresponding ray-theory calculations (dotted curves). All values have been averaged over a range of ±2.4° in travel distance and 
smoothed over ±3° in latitude. Error bars in each panel show the size of the latitudinal average of the standard deviation, cTjvSi 
a measure of the scatter between individual measurements, and the standard error of the mean, ?jvs> as defined by Equations 6 
and 7, respectively. Also, measures of the deviation of the measured travel-time diffrences from the their ray-approximation 
values are show by /ijvs f° r the N-S travel-time differences plotted here, and for comparison the corresponding Hew f° r the 
E-W travel-time differences from Figure 4. 



can be computed by integrating along the ray path 
the flow component tangential to the ray. Specif- 
ically, in the first approximation the travel-time 
difference between waves going along a ray path 
and going along the same path in opposite direc- 
tion is given by: 

At = -2 f ^^ds (4) 

J raypath ^0 

where T is the unit vector tangential to the ray 
path (e.g., Kosovichev & Duvall 1997). The calcu- 
lations in this study used ray paths computed from 
the ray-tracing code of Couvidat & Birch (2009). 

We can see from Figure 3 that the results from 
the analysis of the simulated data are very close to 
the expected travel-time differences computed us- 
ing ray theory. This is true especially for the four 
smaller travel distances. The signal-to-noise ratio 



(S/N) is high here since the deviations from the 
ray-approximation calculations are rather small. 
For the two larger distances, however, the noise 
seems to be comparable with the signal. 

4.2. Error estimation 

When working with measurements, whether us- 
ing actual observations or artificial data from nu- 
merical simulations, it is imperative to estimate 
the accuracy of the measurement results. This has 
been an important issue in helioseismology where 
the results are obtained using complicated data 
analysis procedures. Therefore, we present the er- 
ror estimation in detail. Several error estimates 
are shown in Figure 3 and are explained in the fol- 
lowing. Each data point in Figure 3 is the mean 
of a large number of individual measurements - 
these are the travel-time differences from individ- 
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Fig. 4. — Travel-t imc differences between east-going and west-going waves (solid curves) computed from the same data and for 
the same travel distances as in Figure 3. Values have also been averaged over a range of ±2.4° in travel distance and smoothed 
over ± 3° in latitude. Since no flow in east-west direction was prescribed in the simulation, the travel time differences should 
be zero (dotted curves) if there was no measurement noise. Given error bars are analogous to those in Figure 3. 



ual fittings of the cross-correlation function for the 
different measurement tiles, a range of distances, 
and a range of latitudes. From these, we can com- 
pute a statistical error estimate. Let us denote 
the individual measurements by yi with i <E [1, N], 
where N is the number of individual measurements 
averaged for each final data point, y. Here: 



1 N 



N 



(5) 



The scatter in the individual measurements is 
given by the standard deviation, <r NS , defined as: 



\ 



1 N 



(6) 



(Tns is a measure of the deviation of an individ- 
ual measurement from the mean y. However, we 
are not really interested in the error of an individ- 
ual measurement but rather in the accuracy of the 



mean y. Its error can be estimated by 

&NS 



N 



(7) 



if we can assume statistical independence of the 
individual measurements. This estimate, e^vs, is 
often called the standard error of the mean. As 
can be seen from Figure 3, cms is very small 
and is in fact significantly smaller than the dif- 
ference between the measurements and the ray- 
approximation calculations. This indicates that 
there arc either systematic differences between the 
measurements and the ray-approximation calcu- 
lations or that statistical independence is a poor 
assumption. 

An independent error estimation can be derived 
by looking at the deviation of the E-W travel-time 
differences from their expected values. It seems 
reasonable to assume that the error in both di- 
rections is of the same order. Results for the E-W 
travel-time difference are shown in Figure 4. Since 
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no azimuthal flow is prescribed in the simulation, 
E-W travel-time differences should strictly vanish 
for all latitudes without using any approximation, 
such as the ray approximation. However, due to 
the finite length of the measurement, non-zero val- 
ues are found. The deviations can be used as an 
estimate of the error in both E-W and N-S travel- 
time differences. We define: 



Hew 



1 M 

^E^)-w)(fi)r 



(8) 



where y(0j) is now the mean E-W travel-time dif- 
ference as a function colatitude 9j , M is the num- 
ber of different colatitudes, and yo(&j) is the the- 
oretically expected travel-time difference and is 
zero in this case. As is evident from Figures 3 
and 4, Hew seems to be a more appropriate esti- 
mate of the noise in both the N-S and E-W travel- 
time measurements than <jns and (Jew- Also 
shown in these figures are values of the deviation 
of the N-S travel-time difference from their ray- 
approximation values, i.e. Hns- It is computed 
like equation 8 except that yo is replaced by the 
N-S travel-time difference computed from the ray 
approximation. Values of /ins and hew are sim- 
ilar indicating that the systematic error made by 
the ray approximation is probably small compared 
to the statistical variability of the measurements. 
It is somewhat surprising that hew is actually 
slightly larger than /lijvs for the longest travel dis- 
tances considered here. However, the scatter in 
the N-S and E-W travel-time measurements, i.e. 
<jns and (Tew, are very much the same as seen in 
Figures 3 and 4. 

4.3. Signal-to-noise ratio 

Using the amplitude of the ray-approximation 
N-S travel-time difference and an estimates of the 
measurement noise in the form of the deviations 
Hns, w e can compute a signal-to-noise ratio (S/N) 
for our measurements. The results are shown in 
Figure 5 for a range of measurement lengths. As 
one would expect, S/N increases approximately as 
the square root of the measurement time. We can 
use this dependency to roughly estimate the du- 
ration of a time series one would need to measure, 
say, the flow near the bottom of the convection 
zone to a certain accuracy. For a heliographic dis- 
tance of 42° and 72 hours of measurement time 
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Fig. 5. — S/N of the N-S travel-time difference measure- 
ments as a function of the length of the analyzed time se- 
ries for different travel distances. S/N is here defined as 
the ratio between the maximum amplitude over latitudes 
of the ray-approximation N-S travel-time difference and the 
standard deviation, /i/vSi °f the measured travel-time dif- 
ferences from their ray-approximation values. 
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Fig. 6. — Ratio between the amplitudes of the N-S travel- 
time differences from ray approximation computed for the 
return flow alone, At^p, and for the full meridional flow, 
Atjvs, as a function of travel distance. The peak contribu- 
tion from the return flow in this case corresponds to about 
0.9 seconds in travel-time difference. 
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(3 days) we have a S/N of approximately 1.25 ac- 
cording to Figure 5. However, remember that the 
flow velocities in the present simulation have been 
increased from their realistic values to 500 m/s 
at the surface to make the simulation and mea- 
surements feasible. For a more realistic meridional 
flow of, say, 20 m/s, the signal-to-noise ratio would 
be 25 times smaller. Assuming a S/N of 2 is de- 
sired, one would then need (25 x 2/1. 25) 2 x 3 days 
or on the order of 10 years. MDI mcdium-1 mea- 
surements are available for almost continuous 15 
year, so this may be already possible. However, 
note that the signal for such long measurement 
distances still remains dominated by the strong 
poleward flow in the upper layers of the Sun. As 
shown in Figure 6, only a small portion of the sig- 
nal - up to about 20% at its peak - is from the 
return flow which in the model considered here 
starts at a depth of approximately 146 Mm. So, a 
higher signal-to-noise ration may be desired. 

5. Conclusions 

We have simulated the propagation of acoustic 
waves in the full solar interior in the presence of 
a prescribed meridional flow with a deep return 
flow, and we performed time-distance helioseis- 
mology measurements to detect the effects of the 
meridional circulation on the acoustic travel-times 
difference between north- and south-going waves. 
The measurements were done for large travel dis- 
tances between 12 and 42 heliographic degrees cor- 
responding to lower turning points of the acoustic 
waves between 52 and 195 Mm below the pho- 
tosphere, i.e. deep in the solar convection zone 
all the way to the tachoclinc. The flow velocity 
in the model was artificially increased by a signifi- 
cant factor to a value of 500 m/s in order to model 
the flow measurements using relatively short time 
series that can be calculated on currently avail- 
able supercomputer systems. The results show 
that this approach works well without significantly 
changing the physics of wave propagation, as ex- 
pected from theoretical grounds. The results also 
show that it is, in fact, possible to measure the ef- 
fects of a meridional flow in the deeper solar con- 
vection zone by employing a deep-focusing time- 
distance helioseismology technique. Within the 
statistical variability (noise) of the measurements, 
the measured N-S travel-time differences agree 
well with the ray-approximation calculations. For 



distances between 12° and 30° corresponding to 
lower turning depths between 52 and 136 Mm, the 
agreement is in fact excellent, and still good for 
36° (164 Mm depth). Noise starts to dominate for 
the very longest travel distance, however. We esti- 
mate that for realistic values of the meridional flow 
velocity ~ 10 year time-series or longer may be 
needed for adequate S/N. Such data are currently 
available from the SOHO and SDO space observa- 
tories (since 1996), and ground-based GONG net- 
work (since 1996). It should be mentioned that 
the present simulation uses rather simple mod- 
els for the excitation of acoustic waves as well as 
wave damping, and that therefore the noise prop- 
erties of the Sun may not be very accurately rep- 
resented in this numerical model. None-the-less 
does it seem clear than very long helioseismology 
observations are needed in order to detect small 
flows at the base of the convection zone. Still, 
S/N may be increased, for instance, by the use 
of phase-speed filters which we have not explored 
here, or by more spatial averaging. We also leave 
it for future work to develop and perform an inver- 
sion to infer actual flow velocities from the mea- 
sured travel-time difference. It seems, however, 
that the current ray-approximation based travel- 
time inversion techniques are sufficiently accurate. 
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